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Abstract 

Cosmic strings are topological defects thought to have formed early in 
the life of the universe. If such objects exist, a study of their interaction with 
black holes is of interest. The equations of motion of a cosmic string have 
the form of highly non-linear wave equations. General analytic solutions, 
except for motion in certain backgrounds such as flat and shockwave space- 
times, remain unknown; consequently, much of the work must be carried out 
numerically. To do this, an implicit hnite difference scheme was developed 
that involves solving large block tridiagonal systems. This paper discusses 
the numerical method, its validation against analytic and semi-analytic so¬ 
lutions, and preliminary results on the interaction of a cosmic string with 
Schwarzschild black holes. 


1 


1 Introduction 


This paper [| discusses numerical solutions to the equations of motion of a 
cosmic string in black hole spacetimes. Strings can be either closed (a loop) 
or open. Open strings can be of infinite length, or they can be terminated 
by monopoles; strings can also terminate on the horizon of a black hole (due 
to earlier capture). From the numerical point of view, cosmic strings of 
finite length are preferable since they lead directly to finite computational 
domains. The goal is to model the string as a segment terminated by massive 
end points, taking the limit where the masses go to inhnity to simulate an 
inhnite open string. 

A cosmic string is an extended object under tension; motion of the string 
near a black hole represents the resultant of the competing influences of 
tension and gravity. The results presented here concern the deformations 
induced by the gravitational field of Schwarzschild black holes. These results 
will be compared to the behaviour of dust strings, one-dimensional distribu¬ 
tions of test particles that mimic a string with no tension. 


2 Equations of Motion of a Cosmic String 


The equations of motion for a cosmic string with massive end points are 
derived from an action functional. Solutions to these equations describe the 
world surface of the string as functions .^^) of the string parameters 

= 0,1) conventionally shown as (^°,^^) = (t, cr). The derivation builds 
upon two sources, Larsen and Frolov 0 , who obtain equations of motion for 
open strings in arbitrary backgrounds, and Barbashov 0, who obtains equa¬ 
tions of motion for strings with massive end points in Minkowski spacetime. 
The action for a string terminated by massive end points in an arbitrary 
spacetime is taken to be. 
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^Preprint: Alberta Thy 12-97. To appear in Proceedings of Seventh Canadian Confer¬ 
ence on General Relativity and Relativistic Astrophysics (Fall 1997) 
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where /i is the string tension; cTj(r) represent the motion of the end points 
in the 2D parameter space of the string; Hab is the metric of the string pa¬ 
rameter space, with determinant h; Gab = Qfiv ||x fp- is the induced metric 
on the string worldsheet; and A = A(r) is the most general parametrization 
for the end points whose motion is described by (A(r), c’'i(r)). 

The end points are taken to represent the mass of a longer (potentially 
inhnite) string lying outside the region of interest; the motion of these points 
represents a boundary condition for the section of the string worldsheet under 
study. In the limit of inhnite mass, these end points move along geodesics. 

The equations of motion are obtained by taking the and ai 

variations of the action. The details of this derivation are given elsewhere 
[^]. The internal string coordinates are chosen such that r is a time-like 
coordinate, and a a spatial coordinate used to parametrize points along the 
string. The resulting equations of motion for the interior of the string and 
boundary conditions are. 
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which are used as checks on the quality of the numerical solutions to (|^). 


3 Discretization of the Equations of Motion 

The above equations have the form of non-linear wave equations. Ames P| 
describes a discretization due to Von Neumann for the linear wave equation, 
Utt = Uaa, which uses a standard second-order centered difference formula 
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for evaluating the time derivative at the discrete grid point u (aj, tj), and a 
weighted average (weighting parameter A) of centered spatial differences at 
three adjacent time steps, Tj_i, Tj, and Tj+i. This averaging gives rise to a 
9-point implicit scheme. Von Neumann’s discretization must be extended to 
handle a system of equations and the non-linear terms containing Christoffel 
symbols. To simplify the notation, denote = Xjj. The second- 

order time derivative is discretized using a centered difference formula on a 
non-uniform mesh (to accommodate the iterative scheme described below). 
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while the second-order spatial derivative is discretized using the standard 
centered difference, denoted j, as an average of three time levels. 
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The non-linear spatial term, denoted H for compactness, is also dis¬ 
cretized using the averaging method, (^f)V ~ -|- (1 — 2 A) Hi j + 

XHij-i, where the upper index is an iteration index. The contributions 
at the j and j — 1 level can be obtained directly from centered differences, 
evaluating the Christoffel symbols (using analytic expressions) at Xjj and 
Xj j_i. The j + 1 term is linearized through a Taylor expansion, = 

+ -^(x”^-^^) • Ax”^._,_j^ where Ax”^._,_j^ represents the difference between 
the iterated solutions at the j + I time-level, Ax((^._,_j^ = and 

J (x) is the Jacobian of the n-th iteration of The temporal non-linear 

term is expanded using centered differences on a non-uniform temporal mesh 
and rearranged using the A-notation introduced above. This term is lin¬ 
earized directly by discarding terms of order Ax^. 

The various difference expressions are combined and the terms involving 
the iterated solution, AxOj_,_j^, isolated on the right hand side. The system 
reduces to an equation for the components of a block-tridiagonal system. 


j+i . • Ax“,+, + . Ax“^.,,j^, = 


( 6 ) 


where A, B, and C are 4x4 matrices and d is a 4-vector. The block- 
tridiagonal system expresses an implicit relationship between the iterated 
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solution at all points on the spatial grid (of N points) and reduces to the 
matrix expression A x” = where T in an A x A matrix whose elements 
are 4x4 matrices, and x” and d"’ are vectors of length N whose elements are 
4-vectors. The iterated solution for the new time level is found by inverting 
matrix T; the process is repeated until the solution converges (i.e. produces a 
zero d-vector). Failure to converge triggers a time-step reduction mechanism; 
if this mechanism fails to hnd a suitable new step size, the solution is stopped. 
The constraint equations (|^) are invoked periodically to assess the quality of 
the solution by computing an average value (and standard deviation) for the 
constraints over the spatial grid. 

4 Testing and Using the Solver 

The hnite difference expressions were coded and tested against an analytic 
flat space solution (in spherical-polar coordinates so as to test the non-linear 
portions of the code) and a semi-analytic shockwave solution. The former 
test established the accuracy of the method, the latter tested the iterative 
scheme at the ultra-relativistic extreme. 

4.1 Minkowski Spacetime 

An analytic expression is easily obtained for a cosmic string of length L 
moving at constant velocity v in Minkowski spacetime. In Cartesian coor¬ 
dinates, the string is oriented along the z-axis, such that Z(r, a) = Lcr/Tr, 
(—I < a < ^), with motion taking place along the x-axis. The motion is 
described by, 

x(r, a) = (— cosh (/3) t + Tq,— sinh (/3) r -F Aq, Aq, - cr) (7) 

where Xq, and Tq represent arbitrary initial values. The parameter (3 
represents the rapidity of a string moving with velocity v = tanh/3. This 
solution was coded in spherical-polar coordinates to exercise the non-linear 
portions of the code, and is compared to the numeric solution in Figure 1, 
where the absolute relative error averaged over the spatial grid is plotted as 
a function of string proper time, r. The graph covers roughly 10,000 time 
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steps, and the slope confirms that the method is second-order in the time 
step size At. 

This solution (expressed in Schwarzschild coordinates, (t,r,9,(j))) is also 
used as initial data for all numerical work, since black hole spacetimes of 
interest are asymptotically flat. The initial time to is taken to be zero. The 
value of tq sets the initial distance from the black hole and is chosen to make 
the error of the analytic solution small. The value of 0o sets the initial impact 
parameter b {b = ro sin 0o) relative to the black hole. The string is oriented 
symmetrically about the equatorial plane. 


4.2 Shockwave Spacetime 

The shockwave spacetime arises when an extremely relativistic string passes 
near a black hole. In this situation, the string effectively propagates in 
Minkowski spacetime up to the point where it encounters the sharply de- 
hned gravitational potential of the black hole. 

In the limit where string velocity n —*• c, the Schwarzschild and Kerr 
metrics can be reduced to a shockwave form (Hayashi and Samura, p), 

ds^ = —du dv + dy^ -|- dz'^ — 4:pg{p) 6{u) du^ (8) 


where, 

M = f-x;n = f-Fx;p = 7 Mbh] p= \Jy^ + (9) 

and, for the Schwarzschild case. 
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In this metric, the string is propagating in the x-direction and the coor¬ 
dinates y and z are transverse to the motion. Spacetime is everywhere flat, 
except at M = 0, so the flat-space constant velocity solutions are used as a 
starting point. The equations of motion of strings in a shockwave background 
can be expressed as a Fourier series (adapted from Amati and Klimcik, ||I[), 

U (a, r) = T 

o'" 

V{(T,t) = Vq + p" T + i COs(n(T) 
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where p" is the a-component of the string’s 4-momentum, and the Fourier 
coefficients for the V and Z solutions are given by, 



( 12 ) 


where /(X((j, 0)) = —4pg(X(a,0)), and is evaluated numerically. 

It is possible to compare this series expansion with numerically-generated 
solutions for a string moving at ultra-relativistic velocity. The most interest¬ 
ing comparison is made for projections of the worldsheet onto the Y-Z plane. 
Figure 2 shows the numerically generated solution, the shockwave series so¬ 
lution (to 300 terms), and a numerically generated dust string solution where 
an array of test particles with identical initial data as the cosmic string are 
propagated using the geodesic equation for test particles. 

4.3 Close Encounters 

A cosmic string straying too close to a black hole will become trapped; the 
conditions under which this happens, expressed in terms of the 2D parameter 
space of initial velocity v and impact parameter b, can be used to generate 
a plot of the critical impact parameter (a line in b-v space). These results 
are discussed in a separate paper [Q. The previous section discussed the 
excitation of ultra-relativistic strings; here we consider the case of a string 
moving more slowly. The study of geodesics of particles shows that for impact 
parameters near the critical value, particles can execute a partial orbit around 
the black hole before escaping. Figure 3 shows the trajectory of a slow 
string with near-critical impact parameter, displaying a time-sequence for 
the portion of a string in the immediate vicinity of the equatorial plane, 
which shows the string attempting a partial orbit before being pulled back 
into its original direction of motion. 
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5 Discussion 


In this paper, equations of motion for a cosmic string were derived from an 
action that deals with a finite string terminated by massive end points. The 
end points are taken to represent the mass of segments (potentially infinite) 
lying outside the region of interest; the motion of these points represents 
a boundary condition for the section of the string worldsheet under study. 
In the limit of infinite mass, these end points move along geodesics. The 
equations of motion for the interior of the worldsheet have the form of highly 
non-linear wave equations. 

A numerical method was developed based on Von Neumann’s discretiza¬ 
tion of the wave equation. The discretized equations are linearized and give 
rise to a block-tridiagonal system of equations that must be solved iteratively 
to advance the solution in time. 

The validation of this method against analytic solutions in Minkowski 
spacetime and semi-analytic solutions in shockwave spacetimes was discussed, 
as well as a comparison against dust strings. The excitation of a slowly mov¬ 
ing string with impact parameter near the critical parameter for trapping 
was also shown. 

These results illustrate the role of tension in determining the dynamics of 
the string. At ultra-relativistic speeds, a string develops a strong excitation 
in the direction normal to its original motion. The cosmic string and the dust 
string both exhibit loop formation immediately following the encounter, as 
portions of the strings nearest the equatorial plane cross over to the other 
side. Whereas the portions of the dust string undergoing this cross-over 
continue in their new direction, those portions of the cosmic string that have 
undergone a cross-over eventually recross the equatorial plane, giving the 
worldsheet a kink-like profile that propagates outward along the string. At 
lower velocities, with near-critical impact parameters, loop formation is also 
observed, this time in the direction of motion. The cosmic string performs 
a partial orbit, but eventually gets pulled back in the original direction of 
motion, giving rise to the looping structure. The time sequence shown in 
Figure 3 is for a string with initial impact parameter of 2.28 r^; this value 
is less that the minimum impact parameter of ~ 2.6 for particles (which 
occurs at the ultra-relativistic limit). Numerical studies of the trapping of 
cosmic strings has shown that they are much more difficult to trap than 
particles; this will be discussed in greater detail in |Q. 
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Figure 1: Comparison of numerical solution to analytic solution in Minkowski 
spacetime. 
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Figure 2: Comparison of numerical solution (solid line) to shockwave 
(dashed) and dust (dash-dot) solutions. Initial velocity, 0.995c, impact pa¬ 
rameter 4 Tg. 
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Figure 3: Time sequence of string close encounter. Black hole lies at origin 
of coodinate system. Initial velocity, 0.125c, impact parameter 2.28 Vg. 
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